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Abstract 

in 

A detailed comparison of Faddeev and variational wave functions for 3 H, 
O | calculated with realistic nuclear forces, has been made to study the form of 



three-body correlations in few-body nuclei. Three new three-body correlations 



for use in variational wave functions have been identified, which substantially 
reduce the difference with the Faddeev wave function. The difference between 
the variational upper bound and the Faddeev binding energy is reduced by 
half, to typically < 2%. These three-body correlations also produce a signifi- 
cant lowering of the variational binding energy for 4 He and larger nuclei. 
21.40.+d 
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I. INTRODUCTION 



Quantum Monte Carlo methods, adopted for nuclei in the last decade, have proven to be 
very useful in obtaining exact solutions of the nuclear Schrodinger equation for the ground 
and low-energy excited states of up to six nucleons |[J, and it is likely that many states 
of seven and eight nucleons will be exactly calculated with realistic nuclear forces in the 
near future. A variational approximation, \^f v ), to the lowest-energy state of the desired 
spin, parity, and isospin, (J^T), is first obtained by a variational Monte Carlo (VMC) 
calculation ||. The exact eigenstate, \^o), belonging to the lowest-energy eigenvalue, E , 
for that (J 71 ", T) is then projected out using 

|#o> = Jim exp[-(H - Eq)t]\V v ) (1) 

with the Green's function Monte Carlo (GFMC) method ||. The computation yields a 
population of configurations Rj, where R = (r 1; r 2 , ■ • • , r^) specifies the positions of all the 
nucleons and i — 1, N c labels the configurations, distributed with the probability 

P(R,r) oc \*l{R)exp[-{H - Eo)t]* v {R)\ . (2) 

Mixed estimates, defined as: 

[0{T)]M ~ (^\e X p[-(H - E )r]\^) ' (3) 
are then calculated from this population. Obviously [H(t)]m converges to the eigenvalue 
Eq in the limit r — > oc, however, the variance of [H(t)]m depends primarily on that of the 
local variational energy: 

£ " (R) " *1(r)*.cr) • (4) 

Expectation values of operators other than H are generally calculated only up to first order 
in the difference (|^o) — \^v)) from 

<0> = 2 1ta[0(T)]„-£gg>. (5) 
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It is clear that accurate variational wave functions are needed for the success of this approach. 
The variance of E v (TVj decreases and the terms neglected in eq.(5) become smaller as |\E%) — > 

l*o>- 

State-of-the-art variational wave functions B of few-body nuclei have the form: 



l^,) = [i + E^f + E vsn%) 

i<j i<j<k 



(6) 



with the pair wave function \^f p ) given by: 



I*p> = [^II( 1 + c/ o)]I*j> • 



(7) 



i<3 



The S \\ represents a symmetrized product, and the Jastrow wave function, \^fj), is given 

i<j 

by: 



n/. 



cv i], 



(8) 



Here |$) is an antisymmetric product of single-particle wave functions with the desired 
(J 7r , T), the f c {fij) is a two-body central correlation, and the operator C/y is defined as: 



U, 



E 

p=2,6 



II fijk( T ik, r jk) 



Mm 



u p {rij)0 



(9) 



O?- 1,6 = [1, (Ji ■ (Tj, Sij] <g> [1, Ti ■ Tj] . 



(10) 



The fdjij) and C/y correlations are generated by the static parts of the two-nucleon inter- 
action. The spin-orbit correlations are defined as: 



U ij S = Wlsfaj) + UisrijijYi ■ Tj](L ■ S)ij , 



(11) 



and the eight radial functions fd^ij), w P =2,6( r jj), ^(r^-), and «£ sr (r„) are obtained from ap- 
proximate two-body Euler-Lagrange equations with variational parameters 0. The factor: 

/ \ *2 



Jijk - 1 - £l 



ft. 



ea:p[— £ 3 ft 



(12) 



with 



Rijk ^jk ~\~ T"ki 



(13) 



suppresses spin-isospin correlations between nucleons i and j when a third nucleon k comes 
close to either i or j, and £1-3 are variational parameters. Finally, 



with r = br, represent correlations induced by the three-nucleon interaction Vijk(rij, rjk, r^). 
Its form is suggested by perturbation theory, and e and b are variational parameters. 

The exact ground-state energy of 3 H for the Argonne V14 two-nucleon [|J and Urbana 
model VIII three-nucleon p| interactions, calculated using a Faddeev wave function || and 
Monte Carlo integration, is —8.49(1), while the energy for 4 He, calculated by the GFMC B 
method, is —28.3(2) MeV, where the numbers in parentheses denote the statistical error. 
VMC calculations with the wave function of eq.(6) give instead the energies —8.21(2) 
and —27.23(6) MeV that are too high by ~ 3.5%. This difference, reflecting the deficiencies 
in \^ v ), could be because the two-body functions f c , u p , ii£ S , and it£ ST are not optimally 
determined by the approximate equations used to construct them, or the true |^ ) contains 
many-body correlations in addition to ff^ k and U™ 1 . Note that a \^ v ) containing two- and 
three-body correlations can in principle reproduce the exact |\&o) of 3 H. 

We carried out extensive comparisons between the 34-channel Faddeev \^?f) of the Los 
Alamos- Iowa group ||, one of the best available approximations to the exact l^o); an d the 
for 3 H, to test the accuracy of the two-body functions in \^f v ) and study the three-body 
correlations in \^f), an d therefore in |^o)- These comparisons are discussed in Sec. II. 
Slight improvements in the two-body functions did not significantly reduce the difference 
(\^f) — \^v))- However, a few three-body correlation operators can reduce this difference 
substantially. VMC calculations, for the binding energy of the three and four-body systems, 
with wave functions including the new correlations are presented in Sec. III. The error in 
these is approximately half of that in the best previous variational calculations 0. 



rTNI 
ijk 



(14) 
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II. COMPARISON OF THE FADDEEV AND VARIATIONAL WAVE 
FUNCTIONS AND THE NEW THREE-BODY CORRELATIONS 



The new variational wave function we adopt is given by: 

i*«> = [i+ e tf«*+E^+ e fsrnn /WW. 

i<j<k i<j i<j<k i<j<k 



(15) 



and our purpose is to study the new three-body correlations, f^ k and £4^, and also to 
improve the ffj k correlations in |^ p ). We consider the triton wave function expressed in the 
spin-isospin basis 

IT, 



= \si2sm)\t 12 



22' 



(16) 



where su and tyi are respectively the total spin and isospin of the pair of particles (12), 
s is the total spin of the three nucleons, and m is its projection. Allowed values for s are 
| (for S12 = 0, 1) and | (for s 12 = 1). The total isospin of the system is |, its projection 
is | = — |, and t 12 can have the values or 1. Consequently there are 16 components in 
the wave function, corresponding to the 8 spin states |si 2 sm) labeled by a — 1, 8 and the 2 
isospin states \t\2\\) labeled by (3 = 1,2. These are defined as: 



11/3} = 


10^)1/3) , 


12/3} = 


|oii}|/3> , 


13/3} = 


|li>> , 


K/3} = 


|li>> , 


15/3} = 


|l|2>l/3> , 


16/3} = 




17/3} = 


|l|> , 


18/3} = 


H§> ■ 



(17) 



and 



l«i> = l«>|o^>, 

|a2) = |a)|l~). (18) 

The triton wave function at any configuration R can then be written as 

*(R)=$>«i8(R.)|a/3> (19) 

The main components of the wave function are ^31 and ^12 > which are the only nonzero 
components in the Jastrow wave function. Other components are generated mostly by 
tensor and three-body correlations. 

In order to determine the f^ k and Uijk we compared \^f v ) with \^f), and studied the de- 
viations of the former relative to the latter. Using the decomposition of eq. (|l^) for both wave 
functions, we defined a number of "deviation" quantities. For each spatial configuration, 
Rfc, the deviation £n k is given by 

= E IVv/?( R *) - ^(R fe )| 2 • (20) 

a/3 

Separately, we define a deviation £ a p for each spin-isospin component 

1 N 

^/3=^E I ~ ^F,a/3(Rfc) | 2 , (21) 

k=l 

where N is the total number of configurations considered. Finally, we define an average 
deviation £ 

1 N 

£ = ^E^- (22) 

k=l 

The exact \&o(R) satisfies if\I/o(R) = -E'o^'o(R) a t all R, however the \I/f(R) does not 
satisfy this equation in all regions of configuration space, presumably due to the truncation 
in sums over partial waves and to interpolation between mesh points. In some regions 
^^(R)i7^ F (R) can be different from E ^^(R)^ F (R) by as much as 25%. For the study 
of the deviations listed above, a set of 1,000 configurations R& for which \E f | r (R)-f^ x I / F(R) 
is within 5% of ^o^'f(R)^ / f(R) was chosen from a larger set of random configurations 
distributed with the probability ^(R)^ F (R). 



We studied a wide variety of three-body correlations of the form ^(R)O^. For each 0$, 
values of £ (R^) were obtained by minimizing £ Rfe . Simple functions of R were then chosen 
to approximate the extracted values of £(Rfc). In principle, we can choose a set of sixteen 
operators and define the complete three-body correlation as: 



instead of using eq. (15). The £(R) can then be calculated by solving the matrix equation: 



obtained from ^f v — ty F , at each value of R. However, all our attempts led to unacceptable, 
rapidly fluctuating, functions £ (R) presumably dominated by the small differences between 
typ an d the exact ^ an d inappropriate choices of the sixteen operators 0%. In contrast, 
the £ (R) obtained by minimizing the deviation £ Rfc are smooth and useful for some of the 
operators 0%. 

From the many operators 0$ considered only a few produced a significant reduction of 
the deviations; we list below the successful ones. A spatial three-body correlation, 



reduces the probability for particles i, j, and k to be in a line for q\ > 0. The parameters 
of the new three-body correlations of type x are denoted by qf ; thus q\ and q\ are the 
parameters of f^ k . 

A spin-orbit three-body correlation 



^> = (i + ££(R)^)|* P ) 



(23) 




(24) 



fijk = 1 + Qi( r ij ■ r ik)(*ji ■ r jk )(r ki ■ v kj )exp(-q c 2 R ijk ) , 



(25) 



(26) 



eye 



where 



f'(r) = q[°exp(-qi s r 2 ) 



(27) 



v(r) = exp(—q f 3 s r 2 ) 



(28) 
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may be generated by the two-body spin-orbit interaction u(ry)(L • S)y operating on the 
Jastrow correlations f(r ik ) and f(rj k ). It is expressed as a function of v and / to underline 
its motivation, however v differs significantly from the bare L ■ S potential and both v and 
/ are determined by minimizing the S^ k . 
An isospin three-body correlation 



U ijk = Y,(\ R ijk ~ Tij) [q[exp [-q T 2 (X ljk - g^) 2 ] - 2g[exp(-2g£X 2 fe )} r< • 73 , (29) 

eye " 

with 

Ayfc = [1 + ?X(r»j ■ fifc)(fji • r jk )(r k i • hj)) Rijk ■ (30) 

enhances the probability for = when nucleons % and j are far from k. 

An example of a three-body correlation that we tried and found to be of marginal utility 
is given by: 

U ijk = J2 i(T i' ( T ij X T ik)Uj ■ hk [9ij9ik] m ■ {Tj x r fc ) , (31) 
eye 



where 



g(r) = qfexp(-qfr 2 ) . (32) 

Such a correlation may be generated by tensor interactions between pairs ij and ik. This 
correlation was omitted in our final energy calculations. 



Finally, an improved parametrization of ff jk was found, and eq.([L2|) is replaced by 



fijk = 1 - <?i(l - fife • r jk )exp(-q%R ijk ) . (33) 

Initial values for the various parameters were obtained by minimizing £^ k for the test set 
of 1,000 Faddeev configurations. The final values of the parameters, which are somewhat 
different, are given in Table I; they were obtained by minimizing the binding energy of 3 H 
with the Argonne f 14 and Urbana model VIII interactions for samples of 10,000 variational 
configurations. 
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To illustrate the procedure used to determine the form of the new correlation functions, 
we show in Fig. 1 the data points for the three-body central correlation /^ fe (R) obtained 
by minimizing S^ k at each value of R. These are for isosceles configurations ri3 = r23 with 
Rijk = 6 fm, plotted against cosd, where 9 is the angle formed by ri3 and r23. The functional 
form of eq. fl2"5]) , was chosen because it can fit these numerical values with appropriate 
parameters q\ and q%. The curve in Fig. 1 shows /?- fe (R) with the slightly different values of 
q\ and q\ that minimize the energy; the data points and curve are quite similar. 

The improvement achieved in the main components of ^ v can also be seen in Fig. 2 
where we plot the difference ipF,3i — ip v ,3i for equilateral configurations as a function of R^k- 
The ipF,3i x 10~ 2 is also shown in Fig. 2 for comparison. In this calculation the particles are 
placed in the xz-plane, where ip a p are real, and for equilateral configurations ^31 — — ^12 
holds exactly for antisymmetric The difference ^,31 — i>v,3\ is ^ 1% of ipF,3i for these 
configurations. 

The average values of \ipF,a/3\ 2 , £ a p [eq.(|2l|)1, ^^^1/^ and the total deviation £ [eq.(p^)j 
for the test set of 1,000 Faddeev configurations are listed in Table II for the new and old 
ty v . They are normalized so that ^ F ^p has unit average value for these configurations; 
note that this does not correspond to the standard normalization = (^l^u) = 1, 

however, it clearly exhibits the relative contributions and errors. 

The main components, a/3 = 31 and 12, have a rather small relative error, defined as 
£ap/\^F,a/3 I 2 , of ~ 6 x 10~ 5 in the new \E%. Next in importance is the group of s — | states 
(a = 5 to 8) which have total orbital angular momentum L to t = 2 and account for ~ 12% of 
The largest of these have af3 = 61 and 72, and a relative error of ~ 3 x 10 -4 . The 
states aj3 = 71 and 51 have smaller contributions, and a larger relative error of ~ 10~ 2 . The 
next group has s — |, m — | where the pair of nucleons 1 and 2 is in an odd parity state. 
These states account for ~ 0.1% of \I/p\l/_p and have the largest relative error of ~ 4 x 10~ 2 . 
The last group has s = | and m = — |, for an L tot = 1. It also accounts for ~ 0.1% of 
^p^Fi and these components of the new ty v have an error of ~ 2 x 10~ 2 . 

The £ a p of the new are smaller than those of the old ty v for all values of a(3 other 
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than 62. The total £ of 0.00062 is reduced to 0.00047 by the new f^ k . It is further reduced 
to 0.00035 by the inclusion of Uf? k . The subsequent reduction to £ = 0.00028 is due to U[j k . 
The Ufj k correlations and the changes in the f[j k correlations have very small (< 10~ 5 ) effect 
on the total £ . 



III. BINDING ENERGY RESULTS FOR THE THREE- AND FOUR-BODY 

SYSTEMS; CONCLUSIONS 

As mentioned in the previous section, the final values of the parameters of the new three- 
body correlations were determined by minimizing the energy. In Table III, we present the 
results for the binding energy of 3 H, for the Faddeev wave function and various variational 
wave functions moving from the old to the new. A set of 50,000 configurations sampled from 
(^l^i?) was used for these calculations. The first column contains the values of the binding 
energy and the second the difference Ep — E v . We can clearly see that this difference is 
reduced approximately by a factor of two, when one goes from the old to the new and 
that the f^ k and Uf? k give almost 80% of the improvement. 

The last column of Table III gives (A|A), where 

|A> = |* F > - |tf„> , (34) 

and (^fI^f) = (^f v \^f v ) = 1. The |A) gives the admixture of excited states in \^ v ). The 
mean energy of these excited states is given by: 

TP 

K ~ k nrr ~ 150 MeV (35) 

for both old and new For this reason, the GFMC energy: 

rpf (<K\Hexp\-(H - E a ) T ]\*„) 

|£(T)1 = (^pRiZ-BoMI*,) ' (36) 
of 3 H (and also 4 He) converges to the true Eq by a rather small r ~ 0.02 MeV -1 . In contrast 
the ty v of six-body nuclei seem to have admixtures of states with much smaller E x , and their 
E(t) does not converge so rapidly to E as r increases [|[. 
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These new correlations, obtained by comparing the \I/ f and for the Argonne V\& and 
Urbana VIII Hamiltonian, seem to have more general applicability. Results of variational 
calculations for 3 H and 4 He, with different interaction models, using either the old (Ref. 
[2]) or the new ^> v with the same three-body correlation parameters given in Table I are listed 
in Tables IV and V. (The difference between the old Argonne and Urbana VIII value 
given in Table IV and the old ffj k given in Table III is of marginal statistical significance; it 
can be attributed to slight changes to the the \^f p ) in the present work compared to Ref. [2], 
and the different random walks followed in the two calculations.) A significant improvement 
is seen in each case, with ~ 30-50% of the energy difference with exact calculations, where 
available, recovered. In principle, the parameters of the new correlations should be optimized 
separately for each case, but this has not been done in the present work. The Faddeev 
energies shown in Table IV have been calculated from the Los Alamos-Iowa using Monte 
Carlo integration, and thus have a quoted statistical error. These energies are consistent 
with those published by the Los Alamos- Iowa group ||. For 4 He, the current VMC results 
are fairly competitive with other recently developed methods, Faddeev- Yakubovsky (FY) ]7J 
and correlated hyperspherical harmonics (CHH) ||, as shown in Table V. However, all these 
methods fall short of the exact GFMC calculations in 4 He. 

If the wave function were an exact eigenstate of H, then it should satisfy the eigen- 
value equation: 

H\V) = E\V) (37) 

and this equation should be satisfied for any space configuration, meaning that any local 
energy 

^ (R) " *t(R)*(R) ^ 

should be equal to the eigenvalue E. For the three-body system we considered sets of 
different space configurations, each set being composed of configurations having the same 
Rijk (the perimeter of the triangle formed by the three particles) , and computed the average 
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of the energies returned by eq. fl38|) . These average values, E, are then a function of Rijk 
and can be plotted against it. In Figs. 3, 4, and 5 we show the results respectively for the 
Faddeev wave function, and the old and new variational wave function. These calculations 
were done with the Argonne v u two-body and Urbana model VIII three-body interactions. 

From Fig. 3 we can clearly see that, even for the Faddeev wave function, there are 
regions, in particular for small and large perimeters, where the average values differ from 
the energy expectation value. Moreover, the variance is finite at all R^k, whereas it would be 
zero for an exact eigenfunction. The deviations of E(Rijk) and the variance of the Faddeev 
local energy are perhaps due to the errors introduced by numerical interpolation and the 
application of the permutation operators, required to go from the Faddeev amplitude to the 
full wave function, as explained in the previous section. The application of these permutation 
operators in the Los Alamos-Iowa Faddeev solution is implemented by means of partial wave 
series truncated at 34 channels. This problem led us to consider, in the procedure to search 
for the three-body correlations, only those regions in configuration space where the Faddeev 
wave function gives an E(R) close to the exact energy. A comparison of the old and new 
variational wave function results (Fig. 4 to Fig. 5) shows that the new correlations clearly 
tend to push the average values in the right direction and to reduce the variance. 

In variational calculations of six-body nuclei fl| the new three-body correlations continue 
to lower the variational upper bound. One simple generalization that can be made is to let 
the correlation have different strengths depending on whether none, one, or two of the 
nucleons i, j, and k are in the p-shell. At present it appears best to let this correlation 
be zero if one or two are in the p-shell, although an alternate structure might be even 
better in these cases. It is also beneficial to reduce the overall strength of the Uf? k and U[j k 
correlations. With these steps, the new three-body correlations lower the energy by ~ 0.25 
MeV in 6 Li, and also reduce the variance. 

In conclusion we have reported significant progress in the construction of accurate varia- 
tional wave functions for the nuclear few-body problem. We believe that the new three-body 
correlations will lead to clear improvement in GFMC wave functions for light nuclei, as well 
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as in variational wave functions of heavier systems. 
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FIGURES 

FIG. 1. The three-body central correlation, (f^ k — 1) for isosceles configurations with 
Rijk = 6/m, as a function of cos8, where 6 is the angle formed by fi3 and r23- The dots are 
the numerical values obtained by minimizing Sr and the full line shows the analytical function 
with parameters that minimize the energy. 

FIG. 2. The difference ipF,3i — VV3i> f° r equilateral configurations, as a function of Rijk- The 
dots and the diamonds correspond respectively to the old and new variational wave functions. The 
full line shows ipF,3i(R-ijk) x 10~ 2 . 

FIG. 3. Local energy E(R) of 3 H for 50,000 configurations of the Faddeev wave function, 
sampled with probability \^f\ 2 and binned according to the perimeter variable Rijk', the error bar 
denotes the variance in each bin. The horizontal line shows the average expectation value and the 
curve at the bottom shows the relative distribution of the samples as a function of R^k- 

FIG. 4. Local energy E(R) of 3 H for 50,000 configurations of the old variational wave function; 
notation as in Fig. 3. 

FIG. 5. Local energy E(R) of 3 H for 50,000 configurations of the new variational wave function; 
notation as in Fig. 3. 
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TABLES 



TABLE I. Values of parameters in three-body correlation functions. 



r, 

Qi 


a in 

(J. 19 


1m 


c 
12 


r\ no 


1m 


P 

Qi 


0.20 




j? 

Q2 


0.0b 


tm 


„ls 

Qi 


-O.tz 


tm 




0.12 


fm~ 2 


4 s 


0.85 


fm~ 2 


Ql 


-0.014 


fm _1 


Ql 


0.016 


fm~ 2 


Ql 


1.2 


fm 


Ql 


0.25 
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TABLE II. Average values of the wave functions and deviations. 



a/3 s 12 tl2 S HI \^F,a/3\ 2 £ a (3,new £a/3,old 



31 


1 





1 

2 


1 

2 


4.4 x ltr 1 


2.7 x 10~ 5 


1.2 x 10~ 4 


12 
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1 
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1 

2 


4.4 x ltr 1 


2.7 x 10" 5 


1.1 x 10~ 4 
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3 
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1 
2 


4.1 x 1(T 2 


1.6 x 10~ 5 


1.8 x 10~ 5 


72 


1 


1 


3 
2 


1 
2 


3.0 x 1(T 2 
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2.2 x 10~ 5 


81 


1 





3 
2 


3 
2 


1.3 x 1(T 2 
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52 
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1 


3 
2 
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2 


i.2 x io~ 2 
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62 
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2 


1.2 x 1(T 2 
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71 
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51 
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2 


2.0 x 10~ 3 


1.7 x 10~ 5 
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4.6 x 10~ 4 


1.7 x 10~ 5 
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1 
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4.6 x 10~ 4 


1.7 x 10~ 5 
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42 
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1 
2 
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2 


3.9 x 10~ 4 
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21 








1 
2 
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3.8 x 10~ 4 


5.4 x 10~ 6 


8.7 x 10~ 6 
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2 


1 
2 


2.5 x 10~ 4 


7.4 x 10~ 6 


2.5 x 10~ 5 


22 
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2 


2.1 x 10~ 4 


9.2 x 10~ 6 


2.2 x 10~ 5 


sum 










1.0 


2.8 x 10~ 4 
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TABLE III. Comparison between binding energy values calculated with \^f) and the various 
\^ v ) described in the text. 



Wave function 



E 



E„ — Ef 



(A|A) 



Faddeev 



new f? jk + f? jk + U% k + UT. k 
new fijk + fijk + U ?jk 



new ffjk + fijk 
new f p 



old/f 



ijk 
) 

ijk 



.479 (12) 
.354 (16) 
.346 (16) 
.305 (15) 
.260 (20) 
.253 (19) 



.125 (19) 
.133 (19) 
.175 (19) 
.219 (23) 
.227 (22) 



.00079 
.00099 
.00111 
.00139 
.00146 



TABLE IV. Binding energy in MeV for 3 H with old and new correlations. 









Argonne v 14 


Argonne vi^ 


Argonne V14 


Hamiltonian 


Reid ^8 


Argonne V14 


+ Tucson 


+ Urbana VII 


+ Urbana VIII 


old a 


7.31(2) 


7.45(1) 


8.80(3) 


8.79(1) 


8.21(2) 


new 


7.44(3) 


7.53(2) 


9.05(2) 


8.95(1) 


8.35(1) 


Faddeev 


7.59(2) 


7.70(1) 


9.33(2) 


9.05(1) 


8.49(1) 



■Ref. § 
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TABLE V. Binding energy in MeV for 4 He with old and new correlations. 



Argonne v\4 Argonne v%4 Argonne V14 



Hamiltonian 


Reid i>8 


Argonne V14 


+ Tucson 


+ Urbana VII 


+ Urbana VIII 


old a 


23.62(6) 


23.54(4) 


30.64(9) 


30.51(4) 


27.23(6) 


new 


24.01(8) 


23.80(6) 


31.69(9) 


30.79(5) 


27.63(5) 


FY b 




23.90 








CHH C 




23.93 






27.48 


GFMC d 


24.55(13) 


24.2(2) 






28.3(2) 



a Ref. 
b Ref. § 
c Ref. @ 
d Ref. §] 
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Fig. 1 (Arriaga, Pandharipande, & Wiringa) 



0.02 




cos 



Fig. 2 (Arriaga, Pandharipande, & Wiringa) 



0.024 



0.020 



0.016 



0.012 



0.008 



0.004 



0.000 




